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Abstract 

X 

We introduce an efficient configuration space technique which allows one to com- 
pute a class of Feynman diagrams which generalize the scalar sunset topology to 
any number of massive internal lines. General tensor vertex structures and modifi- 
cations of the propagators due to particle emission with vanishing momenta can be 
included with only a little change of the basic technique described for the scalar case. 
We discuss applications to the computation of ra-body phase space in D-dimensional 
space-time. Substantial simplifications occur for odd space-time dimensions where 
the final results can be expressed in closed form through elementary functions. We 
present explicit analytical formulas for three-dimensional space-time. 



1 Introduction 



Up to now one has not been able to observe any contradictions to the predictions of the 
Standard Model of particle interactions. Possible deviations from the Standard Model or 
revelations of new physics are expected to be quite small at the energies of present accel- 
erators. Future experiments focus on tests of the Standard Model with an unprecedented 
precision |], p|. The present accuracy of experimental data already demands new levels of 
accuracy in the theoretical predictions of perturbation theory J3[. This requirement leads 
to the necessity to compute multiloop Feynman diagrams beyond the one-loop level (as 
a review, see e.g. [|]]). Within the Standard Model the diagrams of perturbation theory 
may contain a multitude of internal and external lines with different particles and different 
masses that makes the task of evaluation of such diagrams rather complicated already at 
the two-loop level. The computation of different subsets of diagrams in different regimes 
of their masses and their external momenta is now an active field of research and often 
requires extensive use of direct numerical methods ||. 

Among the many two-loop topologies the two-loop sunset diagrams with different val- 
ues of internal masses as shown in Fig. 1(a) have been recently studied in some detail 
(see e.g. || [7], |8|, |9|, [10, |ll[] and references therein). In the present note we describe an 
efficient method for computing and investigating a class of diagrams that generalizes the 
sunset topology to any number of internal lines (massive propagators) in arbitrary number 
of space-time dimensions. We call the class of diagrams with this topology water melon 
diagrams. Fig. 1(b) shows a diagram with water melon topology. In our opinion, the 
method presented in the paper completely solves the problem of computing this class of 
diagrams. The method is simple and reduces the multiloop calculation of a water melon 
diagram to a one-dimensional integral which includes only well known special functions in 
the integrand for any values of internal masses. The technique is universal and requires 
only minor technical modifications for additional tensor structure of vertices or internal 
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(a) (b) 




Figure 1: (a) a sunset diagram with three different masses mi, m 2 and m 3 

(b) general topology of the class of water melon diagrams 

(c) one half of a cut n = 6 water melon diagram in the Standard Model representing the 
spectral function for the process q — > qlvpgZH 

(d) an example of a water melon diagram with multiple gluon radiation from internal lines 
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lines (propagators), i.e. tensor particles or/and fermions can be added at no extra cost. 
The method can also handle form factor type processes at small momentum transfer - the 
inclusion of lines of incoming/outgoing particles with vanishing momenta and derivatives 
thereof with respect to their momenta is straightforward and is done within the same cal- 
culational framework. Our final one- dimensional integral representation for water melon 
diagrams is well suited for any kind of asymptotic estimates in masses and/or momentum. 
The principal aim of our paper is to work out a practical tool for computing water melon 
diagrams. In the Euclidean domain the numerical procedures derived from our representa- 
tion are efficient and reliable, i.e. stable against error accumulation. The most interesting 
part of our analysis of water melon diagrams is the construction of the spectral decompo- 
sition of water melon diagrams, i.e. we determine the discontinuity across the physical cut 
in the complex plane of squared momentum. We suggest a novel technique for the direct 
construction of the spectral density of water melon diagrams which is based on an integral 
transform in configuration space. We compare our approach to the more traditional way 
of computing the spectral density where one uses an analytic continuation in momentum 
space. Because the analytic structure of water melon diagrams is completely fixed by the 
dispersion representation, our attention is focussed on the computation of the spectral den- 
sity as the basic quantity important both for applications and the theoretical investigation 
of the diagram. The complete polarization function can then be easily reconstructed from 
the spectral density with the help of dispersion relations. In addition we derive some useful 
formulas for the polarization function in the Euclidean domain and present explicit results 
for some limiting cases where analytical formulas can be found. 
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The paper is organized as follows. In the beginning of Sec. 2 we describe some general 
properties of water melon diagrams and fix our notation. In Sec. 2.1 we introduce the 
configuration space representation of water melon diagrams. In Sec. 2.2 we discuss the 
ultraviolet (UV) divergence structure of water melon diagrams and present the way to 
regularize the UV divergences by subtraction. In Sec. 2.3 some previously known results are 
reproduced using our configuration space techniques. Sec. 2.4 contains explicit examples in 
odd-dimensional space-time. In Sec. 2.5 we discuss expansions in masses and/or momenta 
in the Euclidean domain. In Sec. 3 we consider the computation of the spectral density 
of water melon diagrams by analytic continuation in momentum space. Sec. 4 is devoted 
to the direct computation of the spectral density of water melon diagrams without taking 
recourse to Fourier transforms. Sec. 5 gives our conclusions. 



2 The general framework 

Sunset-type diagrams are two-point functions with n internal propagators connecting the 
initial and final vertex. The sunset diagram shown in Fig. 1(a) is the leading order pertur- 
bative correction to the lowest order propagator in 4 -theory, i.e. it is a two-point two-loop 
diagram with three internal lines. The corresponding leading order perturbative correc- 
tion in 3 -theory is a one-loop diagram and can be considered as a degenerate case of the 
prior example. A straightforward generalization of this topology is a correction to the free 
propagator in n+2 -theory that contains n loops and (n + 1) internal lines (see Fig. 1(b)). 
We call them water melon diagrams^. In a general field theory the internal lines may 
have different masses and may carry different Lorenz structures or may contain space-time 
derivatives. An example of the latter situation is a leading quantum correction in higher 

orders of Chiral Perturbation Theory for pseudoscalar mesons where the vertices contain 
1 Alternative names that have been suggested in the literature for this topology are banana diagrams 
or basket ball diagrams. 
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multiple derivatives of the meson fields. In order to accommodate such general structures 
we represent a general water melon diagram as a correlator of two monomials j n (x) of the 
form 

j n (x)=V lll <f> 1 ---V ltn <f> n (1) 

where the fields <p n have masses m n and where T>^ is a derivative with multi-index fi = 
{/ii, . . . , /ifc} standing for = d k /dx fll . . . dx^ k . The water melon diagrams are contained 
in the leading order expression for the polarization function 

U(x) = (T 3n (x) Jn/ (0)) (2) 

which is explicitly given by a product of propagators and/or their derivatives, 

U(x) = D Mm (ar,mi) • ■■V^ l/n (x,m n ). (3) 

Here V^ u (x,m) = V^V u D(x,m) is a derivative of the propagator D(x,m) with respect to 

the coordinate x with a pair of multi-indices {/i, u}. The propagator D(x, m) of a massive 

particle with mass m in .D- dimensional (Euclidean) space-time is given by 

If e %PtiX ^ d D p [mx) x K\{mx) 
D{x,m) = j—^J p2 + m2 = (27r) A+i x 2A ( 4 ) 

where we write D = 2A + 2. K\(z) is a McDonald function (a modified Bessel function of 
the third kind, see e.g. [0). The propagator D(x,m) depends only on the length of the 
space-time vector \x\ = ^x^x^ for which we simply write x. We consider only the case 
n = n' (equal number of lines at the initial and final vertex). We thus exclude tadpole 
configurations (the leaves of this water melon) which add nothing interesting. Taking the 
limit m — >• in Eq. (^) we get an explicit form of the massless propagator 

D(x0) - 1 r ^dPp _ r(A) 

where T(A) is the Euler T-function. Note that propagators of particles with nonzero spin 
in configuration space representation can be obtained from the scalar propagator by dif- 
ferentiation with respect to the space-time point x. This does not change the functional 



x-structure and causes only minor modification of the basic technique. For instance, the 
propagator for the fermion (spin 1/2 particle) is given by 

d 

S(x,m) = (n^to \-m)D(x,m) (6) 

with 7 M is a Dirac matrix. 

The next generalization consists in extending our configuration space description to the 
case of particle radiation from internal lines. The modification of the internal line with 
mass m due to the emission of a particle with momentum q at the vertex V(p, q) reads 

D{p,m) = — ^ — ? -> — ~V{p,q)- 



p2 _|_ m 2 p2 _|_ m 2 ' — g) 2 4- m 2 

-° -rir^ v M-rir^ = ( 7 ) 

pZ _|_ m z pZ _|_ m z pZ _|_ m ^ 

The functional x-space structure of the corresponding internal line is not changed by this 
modification. Therefore such contributions can be obtained either by differentiating the 
original water melon diagram without particle emission with respect to the mass m or 
by direct differentiation of the propagator which leads to a change of the index of the 
corresponding McDonald function. The explicit representation for the modified internal 
line with mass m is given by 

1 r e ipvxU d D p _ 1 /m\ A ^ . . 

(2^7 (p2 + m 2)/,+i - (2 7 r) A + 1 2AT(/i + 1) \7J ^ mx >- t») 

This modified propagator is a product of some power of x and the McDonald function as 
in the standard propagator shown in Eq. (f4|). The only difference is in the index (A — /i) 
of the McDonald function which is inessential for applications. This does not change the 
general functional structure of the representation constructed below. If there is only one- 
particle emission, the form factor type diagrams with any number of massive internal lines 
can be obtained for any value of q in a closed form. One encounters such diagrams when 
analyzing baryon transitions within perturbation theory. The corresponding formulas will 
be discussed elsewhere. Also the change of the mass along the line in Eq. (0) is allowed for 



vanishing incoming/outgoing particle momenta which allows for the possibility to discuss 
processes with radiation of gauge bosons such as b — > 57 transition for baryons within the 
sum rules approach. As described above, such generalizations can easily be accommodated 
in our approach without any change in the basic framework. From now on we therefore 
mostly concentrate on the basic scalar master configuration of the water melon diagram 
which contains all features necessary for our discussion. 

Eq. (§) contains all information about the water melon diagrams and in this sense is the 
final result for the class of diagrams under consideration. Of some particular interest is the 
spectral decomposition of the polarization function U(x) which is connected to the particle 
content of a model. As a particular example one can consider a water melon diagram 
with a vertex connecting various particles of the Standard Model as shown in Fig. 1(c). 
The knowledge of the analytic structure of the propagators D(x,m) entering Eq. (|]) is 
sufficient for determining the analytic structure of the polarization function U(x) itself. 
For applications, however, one may need the Fourier transform of the polarization function 
n(x) given by 

ft(p) = J U(x)e ip ^dx = J (T Jn (x) Jn/ (0))e ip ^dx. (9) 

The tilde "~" on the Fourier transform Tl(p) will be dropped in the following. The compu- 
tation of an explicit expression for the Fourier transform II (p) for three-line (or two-loop) 
water melon diagrams (the genuine sunset diagram) is described in a number of papers in 
the literature. In the standard, or momentum, representation the quantity H(p) is calcu- 
lated from a (n— l)-loop diagram with (n— 1) integrations over the entangled loop momenta 
which makes the computation difficult when the number of internal lines becomes large. 
Our technique consists in computing the integral in Eq. @ directly using the product of 
propagators Eq. (^) with their explicit form given by Eq. @. The key simplifying obser- 
vation is that the rotational invariance of the expression in Eq. (||) allows one to perform 
the angular integrations explicitly. As a result one is left with a one- dimensional integral 
over the radial variable in Eq. (Q). The remaining integrand has a rather simple structure 
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in the form of a product of Bessel functions and powers of x which is a convenient starting 
point for further analytical or numerical processing. 



2.1 Configuration space representation 

In the present paper we focus on the technical simplicity and the practical applicability of 
the configuration space approach to the computation of water melon diagrams as written 
down in Eqs. @, (f|) and (f|). The idea of exploiting x-space techniques for the calculation 
of Feynman diagrams has a long history. Configuration space techniques were success- 
fully used for the evaluation of massless diagrams with quite general topologies in [13] 



and marked a real breakthrough in multiloop computations before the invention of the 
integration-by-parts technique. 



The case of massive diagrams for general topologies was considered in some detail in |H| 
with x-space integration techniques. As it turns out, the x-space technique has not been 
very successful for general many-loop massive diagrams. The angular integrations do not 
decouple and no decisive simplification occurs. It is the special topology of water melon 
diagrams which makes the x-space technique so efficient. Using x-space techniques one can 
completely solve the problem of computing this class of diagrams. 

So, let us taste the water melon. The angular integration in Eq. (|9|) can be explicitly 
done in D-dimensional space-time with the result 

J d D xe ip ^ = 2tt a+1 X J x (px) (10) 

where p = \p\, x = \x\. J\{z) is the usual Bessel function and dPx is the rotationally 
invariant measure on the unit sphere in the D-dimensional (Euclidean) space-time. The 
generalization of Eq. fllPf ) to more complicated integrands with additional tensor structure 
x \*\ . . . x ^k j s straightforward and merely leads to different orders of the Bessel function after 
angular averaging. The corresponding order of the Bessel function can be easily inferred 
from the expansion of the plane wave function exp(ip At x A ') in a series of Gegenbauer poly- 
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nomials C^p^x^ /px). The Gegenbauer polynomials are orthogonal on a D-dimensional 
unit sphere, and the expansion of the plane wave exp(ip /J x At ) reads 

«p(*jvO = r(A) (^)~ A fv(A + i)j x+l ( P x)cHp^/px). (ii) 

v 1 J 1=0 

This formula allows one to single out an irreducible tensorial structure from the angular 
integration in the Fourier integral in Eq. (P). Integration techniques involving Gegenbauer 



polynomials for the computation of massless diagrams are described in detail in [^] where 
many useful relations can be found (see also fL5|). Our final representation of the Fourier 
transform of a water melon diagram is given by the one- dimensional integral 

n(p) = 2tt a+1 A J x (px)D(x, rm) ■ ■ -D(x, m n )x 2X+l dx (12) 

which is a special kind of integral transformation with a Bessel function as a kernel. This 
integral transformation is known as the Hankel transform. The representation given by 
Eq. (P^) is quite universal regardless of whether tensor structures are added or particles 
with vanishing momenta are radiated from any of the internal lines. An example of a water 
melon diagram with internal gluon emission is shown in Fig. 1(d). 

Next we discuss the analytic structure of a water melon diagram in the complex p- 
plane and also the behaviour of H(p) near threshold. From Eq. ([12]) it is clear that H(p) 
is analytic in the strip |Im(p)| < M = J27=i m i- 111 terms of the relativistically invariant 
variable p 2 this means that the function H(p 2 ) is analytic for Re(p 2 ) > — M 2 implying 
that the function H(p) becomes singular at the energy E = M in the Minkowskian region. 
Depending on the number of internal lines in the diagram this singularity is either a pole 
for the most degenerate case of only one single propagator or a cut in the case of several 
propagators. 

2.2 Regularization and subtraction 



The general formula for the Fourier transform ([12]) allows for an explicit numerical compu- 



tation for any momentum p in the Euclidean domain. However, if D > 2 and the number 
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of propagators is sufficiently large, the integral diverges in the ultraviolet (UV) region or, 
equivalently, at small x. Note that for D = 2 there are only logarithmic singularities at 
the origin (A = for the propagator Eq. (£D) which makes this case technically simpler. 
Also for D = 2 the strength of the singularity does not increase with the number of inter- 
nal lines as dramatically as for higher dimensions. In the general case the UV divergence 
prevents one from taking the limit D — > Dq where Dq is an integer number of physical 
space-time dimensions. Taken by itself, the representation given by Eq. fll2|) determines 
the dimensionally regularized function for complex space-time dimension D. For the scalar 
propagator in space-time with Dq = 4 the singularity at small x is given by x~ 2 . For ex- 
plicit computations in this paper we normally use dimensional regularization, although our 
particular way of regularization is not always the orthodox one. The structure of the UV 
divergence is very simple for the general water melon diagram. Any water melon diagram 
only has an overall divergence without sub divergences if the i?-operation using normal 
ordering and vanishing tadpoles (see e.g. fl6|| ) is properly defined which means that the 
water melon is stripped off its leaves. Thus the renormalization of water melon diagrams 
in the representation given by Eq. flT2"|) is simple and well suited for a numerical treatment 
which is important for practical applications. 

Besides dimensional regularization the momentum space subtraction technique is some- 
times used for dealing with UV divergences. Its implementation is very simple in the rep- 
resentation given by Eq. flT2|). Subtractions at the origin p = are always possible as long 
as one has massive internal lines which prevent the appearance of infrared singularities. If 
this is the case, the subtraction amounts to expanding the function 



(which is a kernel or weight function of the integral transformation in Eq. (|l~2] )) in a Taylor 
series around p = in terms of a polynomial series in p 2 . The order iV subtraction is 




(13) 
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achieved by writing 



px x A 



y ) J\(p%) 



px\- x " (-l) fc fpx\ 2k 



and by keeping N terms in the expansion on the right hand side. Substituting expansion 
Eq. ( |i~4"l) into expression (|l|) leads to a momentum subtracted polarization function 

k 



N 2k / J \ * 

n M o m ( P ) = n( P )-£^— J n( P )| p2 



(15) 



which is finite if the number of subtractions is sufficiently high. The function H(p) itself 
is divergent as well as any derivative on the right hand side of Eq. flT5|) and requires reg- 
ularization. However, the difference or the quantity n Mom (p) is finite and independent of 
any regularization used to give a meaning to the individual terms in Eq. (JT5j). Note that 
the expansion (|T4]) is a polynomial in p 2 in accordance with the general structure of the 
i?-operation |fL6| . The number N of necessary subtractions is determined by the divergence 
index of the diagram and can be found according to the standard rules Hl6| . The subtrac- 



tion at the origin p = is allowed if there is at least one massive line in the diagram and an 
arbitrary number of massless lines. If there are no massive internal lines, the corresponding 
diagram can easily be calculated analytically and the problem of subtraction is trivial. Af- 



ter having performed the requisite subtraction one can take the limit D — > Dq in Eq. ( |12D 
where Dq is an integer. The diagram as a whole becomes finite after the subtraction. In 
order to obtain a deeper insight and for reasons of technical convenience it is useful to 
give a meaning to every individual term in the expansion in Eq. ( |T2D after substituting the 
difference given by Eq. ([14]) which will then finally lead to Eq. ( |15"D . To make the individual 
terms meaningful one has to introduce an intermediate regularization. This intermediate 
regularization can in principle be the standard dimensional regularization. However, the 
regularization in this particular case can be also achieved by adding a factor x 2uJ to the 
integration measure which, in practice, turns into (/^x) 2uJ in order to keep the overall di- 
mensionality of the diagram correct, where [i is an arbitrary mass parameter. Thus, all 
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propagators are taken in the form of an integer number of dimensions (depending on the 
choice of the space-time) but the integration measure is modified to provide regularization. 
We refer to this procedure of auxiliary regularization as an unorthodox dimensional regu- 
larization fl7| . Note that similar modifications of dimensional regularization are known in 
other applications. For instance, in some supersymmetric theories one has to keep the four 
dimensional structure of tensor fields to preserve the Ward identities for the regularized 
theory. The corresponding modification of the standard dimensional regularization is then 
called dimensional regularization by dimension reduction. 

Let us now demonstrate that finite quantities (like the momentum subtracted polar- 
ization operator in Eq. (|l"5D) are independent of the intermediate regularization scheme 
used. Consider the simplest case of one massive and one massless line. The corresponding 
diagram requires only one subtraction according to the standard power counting by noting 
that its divergence index is equal to zero. Within ordinary dimensional regularization the 
polarization function reads 



where 2-^1 (« ; b; c; z) is a hypergeometric function and A = 1 — e. Within the unorthodox 
dimensional regularization we find 



It is straightforward to see that in the limits u — > and e — > both U D (p) — IT D (0) and 
n^o, (p) — n 4ja; (0) are finite and equal to each other in accordance with Eq. (|i5|) and general 
statements of the i?-operation. 

The unorthodox dimensional regularization scheme is only an intermediate step in the 
calculation of finite quantities. However, in some cases it may lead to dramatic simplifica- 




(16) 




(17) 
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tions in the calculation of the relevant integrals. As we shall see later on, the computation 
of the limit D — > D = 3 (and of any odd- dimensional space-time) can be explicitly done 
for any water melon diagram within unorthodox dimensional regularization because of 
the simple form of the propagators and the weight function: they contain the exponential 
functions instead of Bessel functions. The use of ordinary dimensional regularization would 
require the full computation for the noninteger space-time dimension D first which is not 
possible for an arbitrary diagram. 

2.3 Test of the technique with known results 

In this subsection we reproduce some known results with our configuration space tech- 
niques. The complexity of the examples increases with the number of massive propagators 
occuring in the water melon diagram. The completely massless case is quite trivial and 
will not be discussed futher. Water melon diagrams with one massive line and an arbitrary 
number of massless lines are solved with the formula 

jH x»K u (mx)dx = 2^ l m-»- l T + ^ + - \ Y + ^~ - ) (18) 

which allows one to compute all necessary counterterms. Water melon diagrams contain- 
ing two massive lines and arbitrary number of massless lines are also solved in a closed 
analytical form with the formula 

roo 2 2a ~ 3 

/ x 2a ~ 1 KJmx)Kjmx)dx = T(a + {i)T(a)T(a)T(a - n). (19) 

Jo m za L {2a) 

There exists a generalization of Eq. ([19]) to different masses and indices of the McDon- 
ald functions as well. We do not dwell upon this point here. As a simple example we 
give the result for a three-loop water melon diagram with two massive and two massless 
lines at vanishing external momentum. The analytical expression for the diagram in the 
configuration space representation is 

n,o) - /zwzw, = / (<^«) 2 (2o» 
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While the angular integration in D-dimensional space-time is trivial the problem of residual 
radial integration is solved by Eq. (|T5 ). The result for the integral in Eq. (|20|) is 



n(o) 



rn 
~4 



2\ 3A -! 



i r(A) 2 r(i-A)r(i-2A) 2 r(i-3A) 



(21) 



2 % 3A + 3 T(A + l)r(2-4A) 

This result corresponds to the quantity M\ which is the simplest basis element for the 
computation of massive three-loop diagrams in a general three-loop topology considered 
in 
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Next we turn to a sunset diagram with three massive lines, i.e. the two-loop water melon 
diagram. There is an analytical expression for such a diagram at some special values of 
external momenta computed within dimensional regularization |]ID[ . We reproduce this 
result here. We start with Eq. (|^) setting n — n' — 3. The angular integration just gives 
the volume of D-dimensional sphere. We then arrive at Eq. (|12"D with n = 3. Therefore 
the one-dimensional integral to analyze is 

-A 



U D {p) = 2tt a+1 J™ 0?p\ J x (px)D(x, m 1 )D{x, m 2 )D(x, m 3 )x 2X+1 dx. 



(22) 



To localize the finite part we first use momentum subtraction and separate ITd(p) into its 
finite and infinite (but dimensionally regularized) parts 



n D (p) = n M omG») + n sing (p) 



(23) 



where IlMom {p) is a momentum subtracted polarization function and n sing (p) is a countert- 
erm in dimensional regularization. Only two subtractions are necessary by power counting. 
The explicit expression for the momentum subtracted polarization function is 

-A 



n M om(p) = 27T 



A+l 



px 



J\(px) 



mi)D(i, m 2 )D(i, m 3 )i dx 



2ty 



A+l 



1 



p 2 x 2 1 



Jx{px) T(A + 1) ' 4 T(A + 2) 

2A+1 



xD(x } mi)D(x } m2)D(x,m3)x dx. 



(24) 
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The singular part is given by a first order polynomial in p 2 with coefficients A, B whose 
values depend on the regularization scheme (in this case dimensional regularization is used) 



n sing (p) = A + p 2 B 
2tt a+1 



T(A + 



1 poo 

— / D(x ) mi)D(x ) m2)D(x ) m3)x 2X+l 
1) Jo 



dx 



p 2 / x 2 D(x, mi)D(x, m,2)D(x, m^)x 2X+1 dx. (25) 



o 



4r(A + 2) 

With this representation our strategy applies straightforwardly. In the momentum sub- 
tracted part one can forego the regularization (it is finite by /^-operation) and perform 
the one-dimensional integration numerically for D = 4. The counterterms are then simple 
numbers independent of p. They contain divergent parts (regularized within dimensional 
regularization) and need to be computed only once to recover the function Tl(p) for any p. 

In the particular case of the sunset diagram (three lines) the necessary integrals are 
known analytically in their full form and can be found in integral tables (e.g. Using 



the tables may not always be convenient and we therefore present a simplified approach 
which allows one to deal even with the complicated cases in a simpler manner. Let us specify 
to the particular case p = ni\ + rri2 — (pseudothreshold) where an analytical answer 
exists [p70(1 . For simplicity we choose mi = m 2 = m 3 /2 = m. Then p = and IlMom(O) = 



(this is a regular function at this Euclidean point). In this case the counterterm p 2 B 
vanishes because at finite e the quantity B is finite. We thus only have to consider A. 
On the other hand our considerations are completely general since one only requires more 
terms in the p 2 -expansion for arbitrary p. For the special mass configuration considered 
here the result of (ITU reads 



n^ ( o) = ^ m2 " ler2(1 + £) 



3 81n2 , 2 
-- + 81n 2 2 

e 2 e 



+ 0(e). (26) 



(l_e)(l_2e) 

This result can be extracted from the first line of Eq. ([25]) using the integral tables given 
in [ 19| . However, in the general mass caes the necessary formulas are rather cumbersome. 



Even for the special mass configuration considered here they are not so simple. We therefore 
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discuss a short cut which allows one to obtain results immediately without having to resort 
to integral tables. What we really only need is an expansion in e. Basically we need the 
integral 

[■CO 

I D(x,m)D(x,m )D(x,2m)x 2X+1 dx (27) 
Jo 



which is of the general form 



/ x p K ll (mx)K ll (mx)K fJ/ (2mx)dx 

j 



(28) 



(fj, = A and p = 1 — A in our case). For two McDonald functions in the integrand (without 
the last one in the above equation, for instance) the result is given by Eq. (fL9D . Let us 
reduce the problem at hand to Eq. ( JT5| ) and do our numerical evaluations with functions in 
four-dimensional space-time where no regularization is necessary. To do this we subtract 
the leading singularities at small £ from the last McDonald function in Eq. fl28|) using the 
series expansion near the origin, 



r(A) 



i + 



A 



r(i-A) 
r(i + A) {2, 



2X 



0(£ 4 ,£ 



4 f2+2A\ 



(29) 



After decomposing the whole answer into finite and singular parts as (2tc) 2D A = F + S 
(where the total normalization of [|Hj has been adopted in the definition of F and S) we 
find for the singular part 

(2n) D m 2X r°° 
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r(A) 



(2n) D m 
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i-A r(i + A) 
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A-2e 



m 



— T + 



81n2 



(l-e)(l-2e) 
The pole contributions coincide with the result in Eq. 
It is corrected by the finite expression 
(2n) D m 2X 



+ 8(2 -2 In 2 -In 2 2) 



+ 0(s). 



(30) 



r(A + i) 



x 2( - 1 - x) - 1 K x (mx)K x (mx) 
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while the finite part is different. 



(31) 



x { (mx) x K x (2mx) — 



mxf rYl-A), , 2A 



dx. 



1-A r(l + A) 

Because this quantity is finite (no strong singularity at small x) one can put D = 4 to get 

F = 167T 4 m 2 / —K 1 (x)K 1 (x)\xK l (2x)--l + x 2 (-l + 2 lE + 2\nx)\ (32) 
Jo x I 2 L J J 

where 7# is the Euler's constant gamma with the numerical value 0.577. . .. The numerical 
integration results in 

F = 16vr 4 m 2 [-0.306853. . . ] = 167r 4 m 2 [-(l - ln2) x 1.0000. . . ]. (33) 



Now one can restore all e dependence in the normalization factors as in Eqs. (|26"D and ( |30"D 
because F is not singular in e and this change of normalization is absorbed in the 0(e) 
symbol. One obtains 

F = ^~ 2£ {1 _ £){ l ^ [-16(1 - In 2) x 1.0000 . . . ] + 0(e). (34) 

Adding F from Eq. ([34]) to S from Eq. (|30| ) one obtains the result in Eq. (|26"D. Of 
course we use the analytical expression (1 — In 2) for illustrative reasons because we know 
the final answer. 

We emphasize that there is nothing new in computing the polarization function related 
to this diagram at any p. Just some finite part appears (from the momentum subtracted 
polarization function). Also one needs another counterterm for non-zero p 2 . Its computa- 
tion is analogous to what has been done here. It is even simpler because the singularity at 
small x is weaker and only one subtraction from the McDonald function is necessary. This 
technique works for computation at any complex p 2 . 

We add a general remark about the computation near a threshold. The value of the 



polarization function at threshold was obtained in |T0[ . However, the polarization function 
is not analytic in the squared external momentum near threshold so it cannot be expanded 
in a regular Taylor series. Some high order derivatives of the polarization function con- 
sidered as a function of the real variable p 2 do not exist at the threshold. While we do 
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Figure 2: Three loop massive water melon diagram with propagator doubling. The cross 
on the lines corresponds to one differentiation with respect to the propagator mass 



get numerical improvement until the finite derivatives are available (because the very sin- 
gularity is suppressed in these orders), the analytic structure of the singularity cannot be 
obtained in the way of a series expansion. 

The examples presented in this subsection so far are well known and have been obtained 
before using techniques differing from ours. While we can numerically compute any water 
melon diagram with any arbitrary number of internal massive lines it is not easy to find 
corresponding analytical expressions for comparison in the literature. Beyond two loops 
there are only few examples in the literature. We have discovered a relevant result in |L8| 
where an efficient technique for the computation of three-loop massive diagrams of a general 
topology has been developed^. Our test case corresponds to the quantity B N (0, 0, 2, 2, 2, 1) 
in [|18[] and is a three-loop massive water melom diagram at vanishing external momentum 
(see Fig. 2). Three of the four propagators of the diagram are differentiated in the mass, 
i.e. the analytical expression for the propagator in momentum space is (p 2 + m 2 )~ 2 . The 
reference result of |18| reads 

^(0,0,2,2,2,1) = -^C(3). (35) 

The representation in configuration space for such a diagram (after rotation to Euclidean 
2 We thank David Broadhurst for an illuminative discussion of this point 
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space-time) has the form 

wm4 = -4 / K (x) 3 K 1 (x)x 2 dx (36) 
Jo 

where the explicit expression for the derivative of the propagator in Eq. (H) with respect 
to the mass has been used. The diagram is finite (that was the reason for differentiation 
of the propagators). Therefore the expression for space-time dimension D = 4 is used. 
Numerical integration of Eq. ( |3"6"D gives 

wm4 = -2.1036 . . . = -^C(3) x [1.0000 . . . ] (37) 



which coincides with Eq. (|35|). 

We conclude that the configuration space method is simple and rather efficient for 
numerical computation of diagrams with water melon topology. When the structure of 
the transcendentality of the result is known for a particular diagram from some other 
considerations (as in the latest example where only £(3) is present) our numerical technique 
can be used in some cases to restore the rational coefficients of these transcendentalities 
which can be thought of as elements of the basis for a class of diagrams. 

2.4 Examples in odd-dimensional space-time 

It is interesting to note that the computation of Eq. flT2"| ) can be performed in closed form 
for any number of internal lines in space-times of odd dimensions. As the simplest example 
we take three-dimensional space-time D — > Do = 3. 

For A = (A, - 2)/2 = 1/2 with D = 3 the propagator in Eq. (g) reads 

D{x, m) - D 3 {x, m) - - — (38) 

while the weight function after the angular integration given by Eq. ( fL3|) at A = Ao = 1/2 
becomes 

px\- 1/2 T , x 2 sin(px) 

Y) Jl ^ px) = v^^x— (39) 
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The explicit result for the n-line water melon diagram is then given by the integral 

no,, = ^ y,^ ^ 



T(2 -n + 2e) 



(M - ip) n ~ 2 - 2e - (M + ip) n - 2 ~ 2t fi 2t (40) 



2ip(47r) n - 1 

where e is used for regularization and M = J2 m i- 

Here we see the advantage of using the unorthodox dimensional regularization for com- 
puting the finite part of water melon diagrams. The essential simplification of the functional 
form of the integrand within unorthodox dimensional regularization in odd-dimensional 
space-time allows one to compute any water melon diagram analytically. 

We consider some particular cases of Eq. (|40|) for different values of n. For n = 1 we 
simply recover the propagator function with the discontinuity 

= DiscI1 ^ = _L ( H (p) - n(p) , ) = s(s- m 2 ) (4i) 

where s is the squared energy, s = —p 2 . It is usual to call this expression the spectral 
density associated with the diagram. For n = 2 the answer for the polarization function 
n(p) is still finite (no regularization is required) and is given by 

n W = ^.„(^±4 (42) 

Smp \M — ip J 

The spectral density, i.e. the discontinuity of Eq. fl42"D across the cut in the energy square 
complex p 2 plane is given by 

p(s) = ^-^0(s - (mi + m 2 ) 2 ), s = -p 2 , s > (43) 

07l\/S 



which is nothing but three-dimensional two-particle phase space. This can be immediately 
checked by direct computation. The cases with n > 2 have more structure and therefore 
are more interesting. For the proper sunset diagram with n = 3, Eq. ( ^0|) leads to 

' 32tt 2 \e ip \M-ip \ ii 2 I I 
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The arbitrary scale fi 2 appears due to our way of regularization. However, the discontinuity 
of the polarization function in Eq. (f4*4|) is independent of /i 2 as it must be. The discontinuity 
is finite and, therefore, independent of the regularization used. Accordingly, Eq. ([44]) has 
the correct spectral density, 

'M-W^-^- (45) 

The general formula for the spectral density for any n > 1 in D = 3 can be extracted from 
Eq. ©. It reads 

p( 8 ) = ( ^~ M) "" 2 , F 9{s - M 2 ). (46) 



We now want to show how the direct subtraction and our unorthodox dimensional 
regularization are related. Taking Eq. ( f40|) for n = 3 with a subtraction at the origin, one 
obtains 



'o \ px J {Att) 2 x 

which is UV-finite even for e = because there is no singularity at the origin. For practical 
computations it is convenient to keep the factor (n 2 x 2 y in the integrand since this factor 
gives a meaning to each of the two terms in the round brackets in Eq. ([47]) separately. 
Then the direct computation gives 



o y px J (4tt) 2 x 

r(-l + 2e) r ,_ . . . r(2e)//i\ 2 ^ 



2ip(4vr) 2 L v ^ ^ r, } r y M j 

1 \M(M + ip\ , /M 2 +p 2 \l 
1 ln hTT— £ + ln ^T7^ • 



32vr 2 [zp \M-ip) \ M 2 

The poles cancel in this expression and the arbitrary scale \i changes to M. This cor- 
responds to a transition from MS-type of renormalization schemes to a momentum sub- 
traction scheme (with subtraction at the origin in this particular case). Since the spectral 
density p(s) is finite, it can be computed using any regularization scheme as can be seen 
from Eqs. QP]) and 
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We mention that in the three-dimensional case the spectral density p(s) can also be 
found for general values of n by the traditional methods since the three-dimensional case 
is sufficiently simple. The use of the convolution equation needed for the evaluation keeps 
one in the same class of functions, i.e. polynomials in the variable y/s divided by yfs p0| . 
The general form of the convolution equation in D-dimensional space-time reads 

$ n (s) = / M s i)M s *)M s > si, s 2 )dsids 2 , k + p = n. (49) 

For the particular case of three-dimensional space-time the kernel Q 2 {p 2 , ml, m> 2 ) ^ s gi ven 
by 

(2tt) 2 $ 2 (p 2 , ml m 2 2 ) = J 5(k 2 - m\)8{{p - k) 2 - m 2 2 )d 3 k (50) 

or, explicitly, 



Ms, si, s 2 ) = t-^o(s - (VST + V^Y). (51) 

Eq. (^) can be seen to be the two-particle phase space in three dimensions (cf. Eq. ([43"D). 
This is a rather simple example. However, our technique retains its efficiency for large 
n. It is also applicable for odd D other than 3, say D = 5. Namely, the propagator in 
five-dimensional space-time reads (Ao = 3/2) 

(mx) 3 ' 2 K 3 / 2 (mx) e~ mx 

D(x,m) -> D 5 (x,m) = . = ——{I + mx) 52 

(2 , n)°i l x 6 8tt z x 6 

which assures that the integration in Eq. ( p!2|) can be performed in terms of elementary 
functions (powers and logarithms) again. 

We list some potential applications of the general results obtained in this subsection 
for odd-dimensional space-time. In three space-time dimensions our results can be used to 
compute phase space integrals for particles in jets where the momentum along the direction 
of the jet is fixed ||21|| . Another application can be found in three-dimensional QCD which 
emerges as the high temperature limit of the ordinary theory of strong interactions for 
the quark-gluon plasma (see e.g. |9|, ^3[ ^4|). Three-dimensional models are also used 



to study the question of dynamical mass generation and the infrared structure of the 
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models of quantum field theory []25], |26|, ^7|. A further theoretical application consists in 
the investigation of properties of baryons in the limit of infinite number of colors iV c — > oo 
where one has to take into account the spin structure of internal lines. Note that particular 
models of different space-time dimensions are very useful because their properties may be 
simpler and may thus allow one to study general features of the underlying field theory. 
For example, in six- dimensional space-time the simplest model of quantum field theory 3 
is asymptotically free and can be used for simulations of some features of QCD. Though 
five-dimensional models are less popular than others, still there are useful applications for 
Yang-Mills theory in five-dimensional space-time where the UV structure of the models 
can be analyzed ||28|| . 

We conclude this subsection by noting that the x-space techniques allow one to compute 
water melon diagrams in closed form in terms of elementary functions as long as one is 
dealing with odd-dimensional space-times. The resulting expressions are rather simple and 
can be directly used for applications. Having the complete formulas at hand, there is no 
need to expand in the parameters of the diagram such as masses or external momentum. 
In even number of space-time dimensions there are no closed form solutions in the general 
case. In this case it is useful to study some limiting configurations. In the next subsection 
we briefly formulate several limits for water melon diagrams in the Euclidean domain that 
can be computed analytically with the help of standard integrals given in textbooks [29 



2.5 Limiting configurations and expansions 
in the Euclidean domain 

What are the advantages of our method, especially for even dimensions? We shall find that, 
compared to existing approaches, our method results in great simplifications in computing 
the polarization function in the Euclidean domain. While the basic representation in 
Eq. flT2]) can always be used for numerical evaluations, some analytical results can be 
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obtained for particular choices of the parameters in the diagram (i.e. masses and external 
momenta). Different regimes can be considered and some cases can be explicitly done in 
closed form. 

The simplest case is the limit of one large mass and all other masses being small. For 
the polarization function in the Euclidean domain this limit is easy to compute because 
of the simplicity of the x-space representation and the high speed of convergence of the 
ensuing numerical procedures. But this special limit can also be done analytically. When 
expanding the propagators in the limit of small masses one encounters powers of x and 
ln(mx). The remaining functions are the weight function (the Bessel function) and the 
propagator of the heavy particle with the large mass Mh which is given by the McDonald 
function. The general structure of the terms in the series expansion that contribute to 
n(p) is given by 

^ x+l (y) A Jx(px)K u ( y M h x)x 2p \n j ( y mx)dx, p,j>0. (53) 

The integrations in (^) can be done in closed form by using the basic integral representa- 
tion 

^ J\ ip x ) K u ( M h x)dx 

p A r((A + n + v + l)/2)r((A + fj,-v + l)/2) 
2 1 -^M^ +M+1 r(A + 1) 
afi ((A + (i + v + l)/2, (\ + (i-u + l)/2; A + 1; -p 2 /M 2 ) (54) 



where 2*1 ( a > b; c; z) is the usual hypergeometric function [p9|| . The corresponding integrals 
with integer powers of logarithms (j > in Eq. (|53"D) can be obtained by differentiation 
with respect to fi. Note that the maximal power of the logarithm is determined by the 
number of light propagators and does not increase with the order of the expansion: any 
light propagator contains only one power of the logarithm as follows from the expansion 
of the McDonald function K v {^) at small £ . 
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The basic formula Eq. ([12]) is also well suited for finding p 2 -derivatives of the polariza- 
tion function H(p). The values of the polarization function and its derivatives at p 2 = 
can be easily obtained. The convenience of a p 2 expansion is demonstrated by making use 
of the basic formula for differentiating the Bessel function, 

] Jxipx) = {-x 2 ) k Jx+k(px). (55) 



d(p 2 ) k V 2 

Note that differentiation results in an expression which has the same functional structure as 
the original function. This is convenient for numerical computations. Note that sufficiently 
high order derivatives become UV-finite (operationally it is clear since the subtraction 
polynomial vanishes after sufficiently high derivatives are taken). This can also be seen 
explicitly from Eq. (|55|) where high powers of x 2 suppress the singularity of the product of 
propagators at small x. 

Mass corrections to the large p 2 behaviour in the Euclidean domain (an expansion in 
mf/p 2 ) are obtained by expanding the massive propagators in terms of masses under 
the integration sign. The final integration is performed by using the formula 



/ X^J X (px)dx = 2V M 77 TV - 

Jo K ' T (\-^+i\ 



__ , (56) 

2 

Note that all these manipulations are straightforward and can be easily implemented in 
a system of symbolic computations. Some care is necessary, though, when poles of the 
r-function are encountered which reflect the appearance of artificial infrared singularities. 
The corresponding framework for dealing with such problems is well known (see e.g. |]30 



3 Analytic continuation in momentum space 

Now we consider an explicit analytic continuation in the complex p 2 plane and locate 
the discontinuity of the polarization function which is nothing but the spectral density 
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p(s). Note that the spectral density represents n-particle phase space for the water melon 
configuration. 

Again we use the basic formula given by Eq. for analytic continuation. The spectral 
density reads 



2tt Jo 



A+l 



M = if M. M0 



A 



-in(X+l) TT 1 [ m i\T S \ i7r(A+l/2) rr(l) ( m i£ 



i=l 
A 



_ e m(A+l) TT I ( m 'V^ ] e -^(A+l/2) /f (2) ( ^ 

The analytic continuation is performed using the result 



(57) 



where H^' 2 \z) are the Hankel functions, H^\z) = (H^\z))* for real z, A [f29[]. We 
have thus obtained an explicit representation for the spectral density in terms of a one- 
dimensional integral representation. Using Eq. ( |57j ) one can easily discuss special mass 
configurations. It is understood that the UV subtraction has been performed in Eq. ([57]) , 
i.e. the kernel (£/2)~ A J A (£) is substituted with the subtracted one [.]jv from Eq. (|T4]) . 

A remark about the numerical evaluation of the integrals in Eq. (|57D is in order. As 
our main purpose is to create a practical tool for evaluation of the water melon class of 



diagrams, we do not insist on an analytical evaluation of the integrals in Eq. (|57|) . This 
issue will be discussed in more detail in Sec. 4. The one-dimensional integral representa- 
tion in Eq. ( |57j ) is simple enough for further processing. The evaluation of the integral in 
Eq. fllTTD is, however, not always straightforward. The integrand contains highly oscillating 
functions that require some care in the numerical treatment. This is to be expected since 
the discontinuity, or the spectral density, is a distribution rather than a smooth function. 
However, because the analytic structure and the asymptotic behaviour of the integrand 
in Eq. ([57]) is completely known, the numerical computation of p(s) can be made reliable 
and fast in domains where p(s) is smooth enough, in particular far from threshold. One 
recipe is to extract the oscillating asymptotics first and then to perform the integration 
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analytically, or to integrate the oscillating asymptotics numerically using integration rou- 
tines that have special options for the treatment of oscillatory integrands. Both ways were 
checked in simple examples with reliable results. The remaining non-oscillating part is a 
slowly changing function which can be integrated numerically without difficulties. With 
this extra care the integration can be easily made safe, reliable and fast even for an average 
personal computer. We mention that we have checked our general numerical procedures 
in three-dimensional space-time (Dq = 3) where exact results are available (see Sec. 2). 

As an example of the efficiency of our technique we take Eq. Q57D to recompute the 
spectral density of the n = 3 water melon diagram for three-dimensional space-time. The 
Hankel function for indices j+1/2 (or for Dq = 2j+l) is a finite combination of powers and 
an exponential which makes possible the explicit computation of the integral in Eq. (|57|). 
In fact, for this case one has 

1 f°°fsm£ A^fM£\ d£ ^-M n ,_ _ 



p( s ) = / _i _ l S in — ^ — 4— = ^ — e(s-M 2 ). (59) 

HK 1 Trio V ^ / VvW (4vr) 2 e 32tt 2 VS 1 ; K } 

This form coincides with the explicit formulas given by Eqs. (f42|) and fl43|) . The gener- 
alization to higher n includes only algebraic manipulations. The necessary integrations 
corresponding to the one in Eq. (|59|) are performed by moving the contour into the com- 
plex £-plane and regularizing the singularity at the origin by an infinitely small shift ±z0. 
Then closing the contour in the upper or lower semi-plane according to the sign of regular- 
ization one finds the integral by computing the residue at the origin. Note that an explicit 
subtraction is kept in Eq. (|59|). 

The results for the spectral density for n = 3 and 4 can be obtained directly by making 
use of traditional techniques as well. One obtains a one-dimensional integral representation 



for n = 3 and a two-dimensional integral representation for n = 4 pOfl . For larger n, 
however, the corresponding technique of convolution includes many-fold integrals and the 
corresponding recursion relation ( |49"1) is not convenient for applications. This has to be 



contrasted with the one- dimensional integral representation in Eq. Q57| ) derived here which 
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allows one to compute the spectral density for the class of water melon diagrams with large 
number of internal lines in any number of space-time dimensions. 



4 Integral transformation in configuration space 

The analytic structure of the correlator U(x) (or the spectral density of the corresponding 
polarization operator) can be determined directly in configuration space without having 
to compute its Fourier transform first. The dispersion representation (or the spectral 
decomposition) of the polarization function in configuration space has the form 

U(x)= p(m 2 )D(x,m)dm 2 (60) 
Jo 

where for this section we switch to a notation where yfs = m. This representation was 

used for sum rules applications in |33| where the spectral density for the two-loop 



sunset diagram was found in two-dimensional space-time [0. With the explicit form of 
the propagator in configuration space given by Eq. (|]), the representation in Eq. ( |60"D turns 
into a particular example of the Hankel transform, namely the .fT-transform |35], |36| . Up 
to inessential factors of x and m, Eq. ( |60"D reduces to the generic form of the A'-transform 
for a conjugate pair of functions / and g, 



oo 



9(y) = / f(x)K u (xy)y/xydx. (61) 



o 



The inverse of this transform is known to be given by 

] rc+ioo 

f(x) = — g(y)Iu(xy)y/xydy (62) 

HI Jc-ioo 

where I v (x) is a modified Bessel function of the first kind and the integration runs along 
a vertical contour in the complex plane to the right of the right-most singularity of the 
function g(y) []36|] . In order to obtain a representation for the spectral density p(m 2 ) of a 
water melon diagram in general D- dimensional space-time one needs to apply the inverse 
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i^-transform to the particular case given by Eq. (|60|). One has 



(27r) fc+ioo 

m x p{m 2 ) = / U(x)x x+1 h(mx)dx. (63) 



The inverse transform given by Eq. (^3|) solves the problem of determining the spectral 
density of water melon diagrams by reducing it to the computation of a one-dimensional 
integral for the general class of water melon diagrams with any number of internal lines 
and different masses. Compared to the general solution given by Eq. (^) the above form 
is simpler. Below we discuss some technicalities concerning the efficient evaluation of the 
contour integral in the representation given by Eq. (^). The analytic structure of the 
correlator in Eq. (|3|) is now explicit using the representation given by Eq. ( |6"3"D and exhibits 
the distribution nature of the spectral density p(m 2 ) as shown below. 

Note that the expression given by Eq. (Q) can have non-integrable singularities at small x 



for a sufficiently large number of propagators when D > 2 |Tj|. Therefore the computation 
of its Fourier transform requires regularization or subtractions |37[]. The spectral density 
itself is finite (the structure of the water melon diagrams is very simple and there are no 
sub divergences when one employs a properly defined i?-operation [|l(|) and thus requires 
no regularization. In the more traditional approach of direct analytic continuation in 
momentum space the explicit representation for the spectral density is given by Eq. fl5T|) 
where one has taken the discontinuity of the Fourier transform across the physical cut for 
p 2 = —m 2 ± iO. This is an alternative representation of the spectral density, and in some 
instances the latter representation can be more convenient for numerical treatment. 

In the following we present some explicit examples of applying the technique of com- 
puting the spectral density of water melon diagrams on the basis of integral transforms in 
configuration space. 
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4.1 One- loop case 

First a remark about the mass degenerate one-loop case is in order. All necessary integrals 
(both for the direct and the inverse if -transform) involve no more than the product of 
three Bessel functions which can be found in a standard collection of formulas for special 



functions (see e.g. p9|). The spectral density in D-dimensional space-time (for two internal 



lines with equal masses mo) can be computed to be 

, 9x (m 2 - 4m 2 l ) A ~ 1/2 . . 

" (m ) = 2^^mX + l/2)rn m > ^ < M > 

This formula is useful since it can be used to test the limiting cases of more general results. 

The corresponding spectral density for the nondegenerate case with two different masses 
nil an d m 2 reads 

(2^ P (m 2 ) = ^ ( A ; V2 m>mi+ m2? (65) 

4m \ 4m 2 / 



where 



n ' = Wm (66) 



is a volume of a unit sphere in ci-dimensional space-time. Note the identity 

(m 2 — m\ — m\) 2 — Am\m\ = m 2 — (mi + m 2 ) 2 m 2 — (mi — m 2 ) 2 (67) 
which immediately allows one to locate the two-particle threshold. 



4.2 Odd-dimensional case 

For odd-dimensional space-time the representation in Eq. (|60|) reduces to the ordinary 
Laplace transformation. To obtain the spectral density (the function f{x) in this partic- 
ular example) one can use Eq. (|6^) . For energies below threshold it is possible to close 
the contour of integration to the right. With the appropriate choice of the constant c as 
specified above, the closed contour integration gives zero due to the absence of singulari- 
ties in the relevant domain of the right semi-plane. By closing the contour of integration 
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to the left and keeping only that part of the function I v {z) which is exponentially falling 
for Re(z) < one can obtain another convenient integral representation for the spectral 
density when the energy is above threshold. The only singularities within the closed con- 
tour are then poles at the origin (in odd- dimensional space-time) and the evaluation of 
the integral can be done by determining the corresponding residues. These are purely 
algebraic manipulations, the simplicity of which also explain the simplicity of the compu- 
tations in odd- dimensional space-time. For a small number of internal lines n the spectral 
density can also be found by using the convolution formulas for the spectral densities of a 
smaller number of particles (see e.g. PUfl ). For large n the computations described in [20] 



become quite cumbersome and the technique suggested in the present paper is much more 
convenient. 

As an example for the odd-dimensional case we present calculations in three dimensions. 
The dispersion representation for three-dimensional space-time has the following form 

/•OO /-00 g~ mI 

H(x)= / p(m 2 )D 3 (x,m)dm 2 = / p(m 2 ) dm 2 (68) 

JO JO 4:71X 

with the three-dimensional scalar propagator 

D3{X > m) ~ (2n)*/*x ~ W (69) 
One can invert Eq. and obtains 

1 rc+ioo 

2mp(m 2 ) = — : / AnxU(x)e mx dx (70) 
2iri Jc-ioo 



which is a special case of Eqs. (p2|) and (p3|) with 



h(z) = J-smh(z) (71) 
2 V 

where one only needs to retain the e z piece in the hyperbolic sine function. The solution 
given by Eq. ([70|) has the appropriate support as a distribution or, equivalently, as an 
inverse Laplace transform. It vanishes for m < M = Yh=i m i since the contour of inte- 
gration can be closed to the right where there are no singularities of the integrand. Recall 
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that for large x with Re(x) > the asymptotic behaviour of the polarization function U(x) 
is governed by the sum of the masses of the propagators and reads 

U(x) ~exp(-xM). (72) 

For m > M one can close the contour to the left in the complex x-plane and then the 
only singularities of U(x) are the poles at the origin of IT(a;) since U(x) is a product of the 
propagators of the form of Eqs. (§) and (^) . The integration in Eq. ( ffOl) then reduces to 
finding the residues of the poles at the origin. Indeed, 

n g— Mx 

U(x)=J[D 3 (x, mi ) = j—^, (73) 

and Eq. ( |7D|) gives an explicit representation of the spectral density through the polarization 
function in a;-space, 

1 rc+ioo 1 /-c+ioo p (m-M)x 

2mp(m 2 ) = / 4irxll(x)e mx dx = . . . / — dx. (74) 

In closing the contour of integration to the left one computes the residue at the origin and 
obtains (n > 1) 

(m — MY l ~ 2 

2mp(m 2 ) = . -. ' -0(m - M). (75) 

Fy ' (4 7 r) n - 1 (n - 2)! y ' y ' 

Eq. ( [75] ) coincides with the expression (flSp. This also explains the simplicity of the struc- 
ture of the spectral density in odd numbers of dimensions of space-time when traditional 
means are used [0. In five-dimensional space-time Eq. (|70|) is applicable almost without 
any change because the propagator now reads 

(mx)V 2 K 3/2 (mx) e~ mx 

D 5 {x,m) = ' = — —(1+mx). 76 

(27r)°/^x d Stt^x 4 

Compared to the three-dimensional case the only additional complication is that the order 
of the pole at the origin is changed and that one now has a linear combination of terms 
instead of the simple monomial in three dimensions. 
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4.3 Even-dimensional case 

For even- dimensional space-time the analytic structure of U(x) in Eq. (^|) is more com- 
plicated. There is a cut along the negative axis in the complex x-plane which prevents a 
straightforward evaluation by simply closing the contour of integration to the left (with 
Re(x) < 0). The discontinuity along the cut is, however, well known and includes only 
Bessel functions that appear in the product of propagators for the polarization function. 
Therefore the representation (|63|) is essentially equivalent to the direct analytic continu- 
ation of the Fourier transform |[7|] but may be more convenient for numerical treatment 
because there is no oscillating integrand in fl63f ). 

In even number of dimensions one is thus dealing with a genuine fT-transform. We 
discuss in some detail the important case of four-dimensional space-time. For D = 4 
(A = 1), Eqs. (D and © give 

U(x) = [ p(m 2 )D 4 (x,m)dm 2 = [ p^m 2 ) ™^™^ dm 2 , (77) 
J J 4n z x z 

and Eq. (^3|) is written as 

] rc+ioo 

2mp(m 2 ) = — A'K 2 x 2 Ii(x)Ii(mx)dx. (78) 

7TI J c—ioo 



All remarks about the behaviour at large x apply here as well. However, the structure of 
singularities is more complicated than in the odd- dimensional case. In addition to the poles 
at the origin there is a cut along the negative axis that renders the computation of the 
spectral density more difficult. The cut arises from the presence of the functions Ki(rriix) 
in the polarization function II (x). Also the asymptotic behaviour of the function I\{z) is 
more complicated than that of Ii/ 2 {z). In particular the extraction of the exponentially 
falling component on the negative real axis is not straightforward. Incidentally, the fall-off 
behaviour of the function I\{z) on the negative real axis can be taken as an example of 
Stokes' phenomenon of asymptotic expansions (see e.g. []12|). While the analytic structure 



of the representation is quite transparent and the integration can be performed along 

34 



a contour in the complex plane, there are some subtleties when one wants to obtain a 



convenient form for numerical treatment analogous to odd- dimensional case [37 . 

After closing the contour to the left (for m > M) using the appropriate part of the 
function I\{z) we obtain 

C c+ioo 



in I x U(x)Ii(mx)dx (79) 

J c—ioo 

r 2 (U(e in r) + n(e _i V)) K x {mr)dr + 2 J r 2 J\{r)K x {mr)dr 

+ / z 2 U(z)(i7iI 1 (mz)+K 1 (mz))dz+ z 2 U(z)(i7rI 1 (mz) - Kximzfjdz 
JC- Jc+ 

for the quantity entering Eq. (|78|). The contours C + and C_ are semi-circles of radius e 
around the origin in the upper and lower complex semi-plane, respectively (see Fig. 3). 
For practical evaluations of U(e ±in r) the following rule for the analytic continuation of the 
McDonald functions is used, 

^i(e ±l7r O = -MO T i7Tlx(0, Z = mr>0. (80) 

Some comments are in order. The polarization function 11(2;) at z = e ±l7T r is propor- 
tional to the product of propagators of the form 

D 4 (e ±in r } mi) ~ — - [^(m^) ± mlx^r)} . (81) 

It is clear from this equation that the leading singular contribution proportional to the 
product of Ki(rriir) cancels in the sum in Eq. (p^). Also the next-to- leading singular term 
disappears because of different signs in the product. Recall that the small £ behaviour of 
the functions Ki(£) and ii(£) is given by 

tfi(0 = ^ + 0tfln0, J 1 (0 = | + O(e 3 )- (82) 

Let us add a few remarks on the final representation Eq. (|79|) which is in a suitable 
form for numerical integration. We have introduced an auxiliary regularization in terms 
of a circle of finite radius e which runs around the origin with its pole-type singularities. 
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The spectral density is independent of e, and the parameter e completely cancels in the 
full expression for the spectral density as given by Eq. ([791). This must be so since the 
spectral density is finite for the class of water melon diagrams. Eq. ( |79|) contains no 
oscillatory integrands (cf. Eq. ([57])), and the integration can safely be done numerically. 
Thus Eq. (|79"1) is a useful alternative representation for the spectral density. In practice 
the integration over the semi-circles is done by expanding the integrand in z for small z 
and keeping only terms singular in e. The expansion requires only a finite number of terms 
and is a purely algebraic operation. Then the singularity in e exactly cancels against those 
of the remaining integrals. This cancellation can also be done analytically leaving well 
defined and smooth integrands for further numerical treatment. 

Even if the full computation of the spectral function in the even-dimensional case 
described here is straightforward it is nevertheless cumbersome. In order to exhibit the 
essential points in this calculation we illustrate it with a simple and instructive example. 
Consider the calculation of the following integral 

r c+ioo i n z 



C — IQO Z 



■e z dz (83) 



which is rather close in structure to the real case. Due to the singularity at the origin one 
has to treat the vicinity of the origin carefully. We proceed by closing the contour to the 
left 

c+ioo l n z r \ nz rco<lm 

e dz = / — —e dz — / — -e dx 



c-ioo Z 2 Jc c Z 2 Ja X 2 



f°° e x dx r In e + id) , ■ 
-2ni ^ + i — -U + ee^e-^c 



x 2 J-tt e 

= -2m jT ^ + i J] + %-« + In e + iA d( f> 

f°° e~ x dx In e 2ni 

= -2tx% \ 1 OH h2vrzlne + 

Je x 2 e e 

„ . / f°° e~ x dx 1 , \ . „. 

= - 2 "{l —-i- la j- (84) 

The combination in the brackets of the last equation remains finite as e — > 0. Let us now 
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consider this limit in more detail. First we split the integration into two parts from e to 1 
and from 1 to infinity, 

r 00 e~ x dx r°° e~ x dx r 1 e~ x dx 

I — =L — + l — (85) 

Then the first integral is just a number which can be found numerically with high precision. 
In the second integral we can expand the exponent in the integrand and find 

/ — =/ F (l-*+ T — ■) = --l + ln« + 5 (l-«) + .... (86) 

The singularity has the correct form and the finite series converges well. If one takes a 
value 0.1 instead of 1 for the splitting point in Eq. (|85| ) the convergence of the finite series 
will be very fast. This procedure would be used for practical integration in a realistic case. 

Now we give an exact answer for our simple example. After two integrations by parts 
we have 

/OO p~ x r]qr /I \ /"OO 

—^f- = e- e (^+mej -J e- x \nxdx. (87) 

The last integral is finite at e = and for our purpose it suffices to compute it in this limit. 
The result is 

poo 

/ e~ x \nxdx = -"iE (88) 
Jo 

where je is again the Euler's constant. All these manipulations can be easily done with a 
symbolic program. For the original integral one finds 

e~ x dx 1 , \ „ . / , /l , \ 1 



277/ { / — lnej = -2-Ki[e~ e [- +\nej Ine 

= 27ri(l-j E ) at e = 0. (89) 



Thus finally 

r c+io ° In z 



c— too Z 



e z dz = 2m(l - 7 B ) . (90) 



This concludes our discussion of how to treat subtractions in this simplified case. The 
generalization to Bessel functions is straightforward (just expand near the origin and get 
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expressions as in Eq. (jSID). Then the form of the subtraction term depends on the number 
n of propagators in a water melon diagram. Writing down an explicit expression for some 
n is routine and we leave it to the interested user. All the required expansions can be 
performed by a symbolic manipulation program. 

Again for n = 3 the results for the spectral density can be obtained directly by tradi- 



tional means through the convolution equation (f49|). It leads to a one-dimensional integral 
representation for n = 3. In this respect the representation in Eq. (|79D is of the same 
level of complexity while the convolution equation is even simpler because it includes only 
an integration over a finite interval. However, for n = 4 the convolution equation leads 



to a two-dimensional integral representation p0|| . F° r larger n the corresponding tech- 
nique of convolution includes many-fold integrals and the corresponding recursion relation 
in Eq. ((49"|) is not very convenient for applications. In the configuration space approach 
all formulas remain the same regardless of the number n of internal lines and thus this 
approach must be considered to be superior to the more traditional momentum space 
approach. 

For two-dimensional space-time the representation analogous to Eq. ( |79"D is simpler 
because there is no power singularity at the origin but only a logarithmic singularity which 
allows one to shrink the contour to a point (take the limit e — > 0). In this case we obtain 

] poo 

p(m 2 ) = - / r(2U(r) - U(e ln r) - U(e- m r))K (mr)dr. (91) 

7T JO 

For D = 2 we present our results for the cases n = 2 and n = 3. For the one-loop case 
n = 2 one has 

1 r°° 

p(m 2 ) = — / rlQ(mir)I (m2r)K (mr)dr (92) 
27T Jo 



which can be integrated explicitly and results in [19 



This series can of course also be directly obtained by expanding Eq. 
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For the case n = 3 one obtains 



P( m2 ) = 7nZ\2 L r(K (mir)Io{m2r)I (m 3 r) + J (mir) K (m 2 r) I (m 3 r) 



o 



+Io(mir)lQ(m2r)K (rri3r))K Q (rrir)dr. (94) 



Again, this result can be obtained by an alternative method, i.e. by making use of the 
convolution equation given in Eq. (fE|). The 5-function &i(m 2 ) = 5(m 2 — m 2 ) is the 
spectral density of the additional internal line with mass m 3 . $2("^ 2 ) and $2(^ 2 , mi, ?n|) 
are the one-loop spectral density given by Eq. fl65|) . In performing the convolution one 
obtains 

1 r(m— mz) 2 rlq, 

p ( m 2) = _L_ r ds i ( 95) 

(2ir) J( mi +m 2 ) 2 yj((si - m\ - ml) 2 - Am\ml)((m 2 - s x — m|) 2 - Asimfj 

The integrand in Eq. fl9~5|) is singular at the end points and a numerical evaluation requires 
some care. Contrary to this, there are no problems when Eq. fl94j) is used. With modern 
computer facilities, the integral representation given in Eq. fl94|) is more suitable for a 
numerical evaluation than the form given by Eq. (f?5|). 

The most direct and efficient way of numerical evaluation of the spectral density is the 
immediate use of Eq. fl63|). After choosing the contour to be a straight vertical line in the 
complex x-plane with some positive c and making change of the variable x — c + iy the 
integration over y is straightforward. For large number of lines n in a water melon diagram 
the convergence of the integral becomes fast because the integrand falls off as |y|~ 3 ( n_1 )/ 2 
at large \y\. One finds several correct decimal figures of the result rather easily using an 
ordinary personal computer. A good check for the quality of numerical results obtained in 
this way is their independence of c. The results must be the same for any particular choice 
of c. 
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4.4 Threshold behaviour 

Finally we return to the threshold behaviour of the spectral density. Using the results of 
the above analysis one arrives at expressions of the general form 

/oo ^ 
Ki(mr)Y\_Zi(mir)dr (96) 

i 

where Z\{miT) is either l\(m,ir) or K\{mir). The convergence at large r > (at the upper 
limit of the integral) is controlled by the factor exp (— (m — M)r) as in Eq. (|74|), and the 
corresponding expansions in the variable m — M in the region m ~ M can be constructed. 
As shown before, the spectral density for odd-dimensional space-time is known exactly. Its 
threshold behaviour is easily extracted and one obtains 

, 9x (m-M)™- 2 , x 

P(m ) ~ — 7 ; 1 j tt 97 

HK ' 2m(4vr)"- 1 (n - 2)! v ; 

at m close to M. 

In the realistic case of four-dimensional space-time the threshold behaviour of the spec- 
tral density can be inferred from Eq. (p3|) (and, also, Eq. (|9~6|)). Substituting the asymptotic 
limit of modified Bessel functions at large arguments we infer from the simple dimensional 
considerations that 

p{m 2 ) ~ (m-M) (3n " 5)/2 (98) 

at m close to M, where n is the number of internal lines of the water melon diagram. 
One obtains p{m 2 ) ~ \Jm — M for n = 2 and p(m 2 ) ~ (m — M) 2 for n = 3 at m close 
to M. The threshold behaviour agrees with the one extracted from the explicit formula 
given by Eq. (^) or with the threshold behaviour derived from the convolution equation 
in Eq. fl49|) . Other space-time dimensions can be analyzed along the same lines. 

For the sake of completeness we briefly comment on the opposite limit for the spectral 
density when m 2 — > oo. In this limit the external energy is much larger than all masses. The 
limiting behaviour of this kinematic configuration can be found by utilizing the massless 
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approximation for the correlator where analytical expressions are available. Subleading 
corrections can be obtained using known techniques ||. 

5 Conclusion 

We have described a novel technique that reduces the computation of any water melon 
diagram to a simple one- dimensional integral with known functions in the integrand. Ad- 
ditional tensor and form factor structures can be easily included without modification of 
the basic formulas. Different regimes of behaviour with respect to mass/momentum ex- 
pansions can be easily analyzed and the results can be found analytically if no more than 
two dimensionful parameters (one mass and momentum or two masses) are treated exactly 
and other parameters are considered to be small and are treated in series expansions. Ex- 
plicit analytical formulas (with even that last integration being performed in closed form) 
are given for the case of any odd number of space-time dimensions. The analytic con- 
tinuation to the Minkowskian region, i.e. to the positive s-axis has been performed and 
the discontinuity across the physical cut is explicitly given. This allows one to compute 
general n-particle phase space for any kind of particles both with different masses and 
Lorentz structures. The threshold behaviour of the spectral density can be easily investi- 
gated based on such a representation. In the even- dimensional case the final integration for 
an arbitrary water melon diagram was not done explicitly since we encountered products 
of Bessel functions which we were unable to integrate in a closed form. This forces one to 
use numerical integrations in the even-dimensional case if numerical results are desired for 
large n. However, the analytic structure of the solution for the diagram is fully determined 
which makes the numerical treatment rather reliable and fast. In this sense the compu- 
tation of water melon diagrams has been converted to the routine procedure of getting 
numbers. Except for some special cases the final one- dimensional integral representation 
was not further reduced to known simpler functions. This is quite familiar from the theory 
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of special functions where many special functions are represented through one-dimensional 
integrals which cannot be further reduced to known functions. In this sense our formulas 
completely solve the problem of computing the class of water melon diagrams. 

Well, the water melon is finished and the sunset is over long ago. So the next story has 
to wait until the next sunrise takes place. 
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